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The issues related to moving elements in space and instruments working in broader wavelength ranges 
lead to a need for robust polarimeters, efficient on a wide spectral domain, and adapted to space condi¬ 
tions. As part of the UVMag consortium, created to develop spectropolarimetric UV facilities in space, 
such as the Arago mission project, we present an innovative concept of static spectropolarimetry. We 
studied a static and polychromatic method for spectropolarimetry, applicable to stellar physics. Instead 
of modulating the polarization information temporally, as usually done in spectropolarimeters, the mod¬ 
ulation is performed in a spatial direction, orthogonal to the spectral one. Thanks to the proportionality 
between phase retardance imposed by a birefringent material and its thickness, birefringent wedges can 
be used to create this spatial modulation. The light is then spectrally cross-dispersed, and a full-Stokes 
determination of the polarization over the whole spectrum can be obtained with a single-shot measure¬ 
ment. The use of Magnesium Fluoride wedges, for example, could lead to a compact, static polarimeter 
working at wavelengths from 0.115 }im up to 7 jim. We present the theory and simulations of this concept, 
as well as laboratory validation and a practical application to Arago. © 2015 Optical Society of America 

OCIS codes: (120.5410) Polarimetry, (260.5430) Polarization, (260.1440) Birefringence, (120.6085) Space instrumentation, 
(350.1260) Astronomical optics, (350.1270) Astronomy and astrophysics. 
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1. INTRODUCTION 
A. Context 

Analysing the polarization state of the light provides crucial 
information in many different areas, for example in remote 
sensing for medical or biophysics issues, military applications, 
atmosphere studies, and of course astronomy [1]. The light 
that we receive from objects in our universe can be polarized 
through several physical processes. Polarized light can be ob¬ 
served from, e.g., the interstellar medium in our galaxy, other 
galaxies, quasars, or supernovae. In stellar physics, the Zeeman 
effect allows us to determine the polarization state of the light 
over the spectrum, which enables the measurement of the mag¬ 
netic field of stars and provides constraints on circumstellar en¬ 
vironments. For observations of the Sun, both the Zeeman and 
Hanle effects are used to determine the local magnetic field in 
and around the Sun. 

In this work, we aim to measure and characterise magnetic 


fields of stars and polarization from the stellar surroundings, 
e.g., the disk or stellar wind. However, it is not yet possible to 
resolve the surface and surroundings of stars in the same way 
as is done for the Sun, with the stellar target apppearing only 
as a point source in observations. Nevertheless, analysing the 
polarization of the light emitted by the star in its photospheric 
absorption spectral lines offers an easy way to measure the 
strength of the longitudinal stellar magnetic field, as it is 
directly proportional to the intensity of the polarized light. The 
presence of a magnetic field manifests as characteristic Zeeman 
signatures in the polarized part of the light, specifically in 
the spectral photospheric lines. Measuring this magnetic 
signature as the star rotates allows us to characterise the full 
magnetic field at the surface of the star and, in particular, its 
polar strength, obliquity with respect to the rotation axis, and 
surface configuration. Similarly, polarization in wind-sensitive 
UV resonance lines, or in emission lines originating from 
disks or choromospheres, allows us to study the circumstellar 
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environment. While some science cases may be tackled by 
measuring circular (Stokes V) polarization only, others require 
linear (Stokes Q and U) polarization. However, it has been 
shown that 3D magnetic mapping of stars and their environ¬ 
ment is more precise when all Stokes parameters are used 
simultaneously [e.g., 2, 3]. In addition, to allow a space mission 
to address various science objectives, it is desirable to provide 
measurements of both circular and linear polarization. 

To characterize and quantify the polarization of light, we 
need to measure its polarization state all over the entire spec¬ 
trum, and with a high spectral resolution. However, as com¬ 
mon detectors are only sensitive to the intensity of the light and 
not to its polarization, we have to build polarimeters to encode 
the polarization state information into the intensity of the light. 
The next section describes the method to determine the polar¬ 
ization states of the light. 


B. Formalisms and tools 

B. 1. Stokes and Mueller formalisms 

To go into more details, the Stokes and Mueller formalisms are 
briefly introduced in this section. They are well documented in 
the literature, for example in the books "Polarized light" from 
W. A. Shurcliff [4] or "Birefringent Thin Films and Polarizing 
Elements" from I. J. Hodgkinson [5]. 

Considering a monochromatic electromagnetic wave prop¬ 
agating in the ^ direction, the transverse components of its 
electric field ^ are defined by: 


Ex = ^x-cos{ojt-kz) 
Ey = sin(a;f -kz + (p) 


with CO its angular frequency, t the time, and k the wave num¬ 
ber related to its wavelength. 

The Stokes vector, giving complete information of the light 
(intensity and polarization state), is then defined using the am¬ 
plitudes ^x and ^y of the electric field and the phase shift cp be¬ 
tween its two components: 


/ \ 

Q 

u 
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The first Stokes parameter I defines the total intensity of the 
light, while the intensity of the polarized part of the light is de¬ 
fined by The ratio of these two values defines 

the degree of polarization of the light p: 


+ 1/2 + y2 

I 


(3) 


The Mueller (4x4) matrix is introduced to link the input and 
output Stokes vectors after an optical element. Adding several 
optical elements in the system corresponds to multiplying the 
Mueller matrices of the elements. Therefore, the input Stokes 
vector we want to determine is linked to the output Stokes vec¬ 
tor we measure with our polarimeter through the total Mueller 
matrix of the system Mtot’ 


^out — ^tot ' 


(4) 


Developing a polarimeter consists in finding the optimal 
combination of optical elements to determine the input Stokes 
parameters from the measured output ones. As detectors are 
only sensitive to the intensity, the goal is that the Mueller ma¬ 
trix transfers the full input Stokes vector into the first Stokes pa¬ 
rameter I of the output Stokes vector. Mathematically speaking, 
we aim at a final equation of the type lout = + ^Q-Qin + 

Xu-Uin + xy.Vin, with Xj, xq, Xjj and xy combinations of some 
elements of the total Mueller matrix, depending on the char¬ 
acteristics of the optical elements. As a consequence, a modu¬ 
lation of the polarization states is needed, because we have a 
single equation with four unknowns, I, Q, U, and V. 

B.2. Modulation of the polarization state 

A typical solution is to modulate the polarization state tempo¬ 
rally, by rotating some optical elements in the system to solve 
the above equation. However, in space missions, rotating or 
more generally moving parts should be avoided if possible in 
the instrument, as they are automatically considered as single 
point failure (i.e., risk items for the mission). Using liquid crys¬ 
tal switching systems could be a solution to modulate tempo¬ 
rally the polarization state without rotating optical parts [6]. 
Unfortunately, liquid crystals are not efficient in wide spectral 
ranges and cannot be used in the far-ultraviolet (FUV) region, 
below 250 nm. In addition, temporal modulation provides suc¬ 
cessive rather than simultaneous measurements of the Stokes 
Q, U and V parameters. Instead of temporal modulation, polar¬ 
ization can be modulated spectrally [7] or spatially using chan¬ 
neled polarimetry or division of amplitude. 

Our goal is to create an efficient static spectropolarimeter 
over a wide spectral range, which provides the full Stokes vec¬ 
tor simultaneously. This means we need to modulate the polar¬ 
ization spatially, over the wavelength range. As we use a spatial 
dimension for the modulation, this method can only be applied 
to point sources, such as stars. This solution imposes a variable 
retardance to the beam in one spatial direction to obtain several 
polarization states at the same time, without moving any part 
of the system. Ideally, the spatial direction in which the polar¬ 
ization will be modulated must be orthogonal to the spectral 
direction, so that we obtain a 2D array containing all the polar¬ 
ization and spectral information of the entrance light. 

Birefringent wedges, generally known as polarization scram¬ 
blers, can be used to create this spatial modulation [8]. Com¬ 
bined with a linear polarizer and a spectrograph, we can create 
a static spectropolarimeter working over a wide spectral range. 
The only limitations of this system are the efficiencies of the ma¬ 
terials and the design of the spectrograph. 

B.3. Material 

Materials such as quartz, sapphire, or calcite are good can¬ 
didates to create a polychromatic spectropolarimeter, because 
of their transparency and birefringent characteristics. How¬ 
ever, for an instrument covering a wavelength range from FUV 
(0.115 ptm) to infrared (IR) (7 f/m), the only available transpar¬ 
ent and birefringent material is Magnesium Fluoride (MgF 2 ). In 
this work, we thus consider MgF 2 wedges. 

The birefringence curve of MgF 2 over the wavelength range 
0.115 to 7 Jim is displayed in Fig. 1. 

The measurements for this plot were obtained by Chan- 
drasekharan in 1968 and 1969 [9, 10] and Dodge in 1984 [11]. 
For the space mission CLASP, aiming at solar spectropolarime- 
try at Lyx, a short part of the birefringence curve has been mea¬ 
sured again in 2013 [12]. The various data were combined to 
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Fig. 1. Birefringence curve of MgF 2 from FUV (115 nm) to IR 
{7iim) [9-12], 


obtain the complete plot of the birefringence of MgF 2 over the 
[0.115:7] Jim spectrum (Fig. 1). 

The birefringence of MgF 2 , An, defined by the difference 
between the extraordinary refractive index Ue and the ordinary 
one Ho, only varies between 0.005 and 0.015 from 122 nm to 
7 Jim. It is transparent down to 115 nm, but its birefringence 
drops off dramatically below 140 nm, goes to zero around 
119.5 nm, and inverts its sign below this wavelength. The 
consequence is that it is very difficult to use MgF 2 to measure 
polarization states between 119 and 120 nm, as the value of 
the birefringence is too low. Nevertheless, MgF 2 offers the 
best opportunity to create a polychromatic spectropolarimeter 
including the FUV domain. 

In the next section of this paper (Sect. 2), we describe the 
conceptual theory, calculations, and simulations of a spectropo¬ 
larimeter using spatial modulation with MgF 2 wedges. After 
applying the method for the design of an instrument for the 
Arago space mission project (Sect. 3), the laboratory implemen¬ 
tation section (Sect. 4) shows the practical results and compari¬ 
son with simulations. 

2. THEORY AND SIMULATIONS 
A. Global principle 

To achieve a spatial modulation of the polarization, we need 
to use a component that gives a variable phase retardance be¬ 
tween the extraordinary and the ordinary part of the light. As 
shown in Eq. 5, the retardance of a plate is directly proportional 
to its thickness. If the thickness d of the plate varies continu¬ 
ously in a certain direction, the retardance O also varies accord- 
ingly. 


^ 2nAn , 
O = —;— • d 


(5) 


We use a birefringent wedge (polarization scrambler) of 
apex angle ^ and birefringence An = Ue — Ug. The wedge is 
illuminated by a collimated beam coming from a point source. 
The link between the retardance and the height x along the ver¬ 
tical direction of the wedge is given in Eq. 6, with x = 0 at the 
top of the wedge: 


^ ,61 

A 

To ultimately create a spectropolarimeter, we then disperse 
the light using, for example, a grating, which causes depen¬ 
dence of the retardance on both the height x and the wavelength 

A. 



Fig. 2. Global principle of the measurement. The light travels 
through two double-wedge birefringent components, a linear 
polarizer, and finally a grating to create the spectrum. A 2D im¬ 
age is obtained with the spectral information in the horizontal 
direction and the polarization modulation in the vertical one. 


Combining a single birefringent wedge with a linear polar¬ 
izer cannot measure the full Stokes vector, but only two of the 
three Stokes parameters in addition to the intensity I (see more 
details in Sparks et al. (2012) [8]). To solve this problem and 
obtain a full Stokes determination, we need to add a second 
birefringent prism with a slope (or apex angle) twice as large 
as the one of the first wedge. In order to make a rectangular 
component, which is easier to handle, and to introduce a sym¬ 
metry around the optical axis, the two wedges can be doubled 
with the same material, but with an orthogonal fast axis angle. 
In this case, on the optical axis, there is no birefringent effect 
and the retardance is zero. This configuration is shown in the 
bottow left part of Fig. 2, and in more detail in Fig. 3. 

As detailed in the next subsection, this configuration of 
wedges followed by a linear polarization analyzer enables us 
to encode the polarization information into the output intensity. 
The light is then dispersed with a grating and focused onto a 
detector. The resulting image contains all the needed informa¬ 
tion: a classical spectrum with, for each wavelength, the full 
Stokes information encoded in the orthogonal direction of the 
spectrum (see Fig. 2). 

It is possible to use a dual-beam analyzer instead of the lin¬ 
ear polarizer. In this case, the signal-to-noise ratio (SNR) in¬ 
creases by a factor of \/2, and some systematic errors (due to il¬ 
lumination inhomogeneity for instance) can be suppressed. Un¬ 
fortunately, such a system would double the size of the required 
detectors, which is a critical issue as shown in Sect. 3. 

The SNR obtained with this concept can be compared to 
other techniques and, in particular, to a more usual polychro¬ 
matic rotating modulator. For this comparison, we consider 
the dual-beam option, as it is a standard technique for spec- 
tropolarimeters with rotating plates. The polychromatic rotat¬ 
ing modulator we consider is composed of 3 quasi-zero order 
MgF 2 plates at 6 different angular positions, optimized to maxi¬ 
mize the extraction efficiencies over the spectral range [13]. For 
the static spectropolarimeter presented here, the flux at each 
wavelength is spread along a column of pixels, while for the ro¬ 
tating system it is divided in several sub-exposures. The ratio 
of the SNR resulting from the two spectropolarimeters is equal 
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Fig. 3. Configuration of the birefringent wedges to obtain full 
Stokes measurement. The first bloc of 2 wedges has a total retar- 
dance (p = (pi—(p 2 arid the second bloc has twice this retardance 
: 2(p. The fast axis angle values (with respect to the horizontal) 
for the four wedges are [^; 0; from the left to the right 

(shown with red dashed-dotted lines). 


to the square root of the ratio between the number of pixels 
used for the concept presented in this article and the number of 
sub-exposures for the rotating technique. 


For the second component, the calculus is the same with 
D 2 = Dn • Dq. The result is directly shown here: 


0 

0 1 

0 0 

^0 0 


0 0 ^ 

0 0 

cos 20 sin 20 

— sin 20 cos 20 y 


( 10 ) 


The global Mueller matrix D of the polarimeter modulator is 
developed in Eq. 11: 


B. Mueller matrix calculation 

Using Mueller matrices applied to this system, this section 
presents in more detail how the polarization is modulated and 
encoded. 

Fig. 3 defines precisely the configuration of the birefringent 
wedges. The fast axis angle of the first wedge of retardance 
Oi is set to a:i = and its apex angle is The fast axis 
angle of the second wedge of retardance O 2 is set to The 
total phase retardance of the first rectangular component is 
then O = — O 2 . Because the two wedges have orthogonal 

fast axes, the retardance cancels for equal thicknesses, and the 
global retardance on the optical axis is thus 0 = 0. 

The second rectangular component is built in the same way, 
but with the fast axis angles of 0° and ^ and wedges with apex 
angle 2^. 

The retardance O is directly calculated using Eq. 6: 


<l>(x, A) = 27r tan ^ • X - 27r tan ^ 

= 271 tan ^ . (2x — d) 

A 


{d — x) 


(7) 


D = D2 • Di 

^1 0 

0 cosO 

0 — sin O sin 20 

^0 — cos 20 sin O 


0 0 ^ 

0 sinO 

cos 20 cos 20 sin 20 

— sin 20 cos O cos 20 y 


( 11 ) 


The final step to determine the equation of the resulting in¬ 
tensity lout after passing through the polarimeter is to compute 
the total Mueller matrix of the polarimeter and to mutliply it 

by the entrance Stokes vector (j Q U The matrix D 

has to be multiplied by the first line of the Mueller matrix of 
the linear polarizer 0.5 • (l cos 26 sin20 O), where 0 is the 

angle of the polarizer with respect to the horizontal. The 2D 
output intensity, depending on the height in the wedge x, the 
wavelength A, and the input Stokes parameters I, Q, U and V, is 
described with Eq. 12: 


hut (^/ A) = 0.5 • [J + Q • (cos O cos 26 — sin O sin20 sin 26) + 
U cos 20 sin 20 + V • (cos O sin 20 sin 20 + sin O cos 20) ] (12) 


The above considerations are made with x = 0 at the top of 
the wedges. Changing the variable to have x = 0 on the optical 
axis (as shown in Fig. 2) leads to a global retardance of: 

0(x,A) = 4/1 tang > ^ (g) 

A 

The Mueller matrix of the first rectangular component Di is 

the product of the Mueller matrices of the first wedge Dstt and 

4 

the second wedge Dn : Di = Dn • Dstt . 


This equation, in combination with Eq. 8, provides the start¬ 
ing point for simulations and to the design of our laboratory 
implementation. 

C. Simulations and inversion 

For the simulations, MgF 2 birefringent wedges, with an apex 
angle ^ = 1.5° are used, as well as a fixed linear polarizer with 
angle 0 = 64°. The value of this angle was chosen so that the 3 
basic 2D patterns corresponding to 100% Q, 100% U and 100% 
V entrance polarization are different. This is not mandatory, but 
makes the demodulation process easier. 
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Fig. 5. Simulation of a 2D image obtained for a synthetic star, 
between 655 and 670 nm. Photospheric lines are observed, but 
the polarization modulation is not seen by eye. 


Fig. 4. Simulation of 2D images obtained by solving Eqs. 8 and 
12 for given polarization states. The values of the constants 
used in this simulation are an apex angle of ^ = 1.5° and a 
linear polarizer angle of 6 = 64°. 


A given entrance polarization state of the light is set: 1=1 and 
Q=l, U=0, V=0 in the first case, Q=0, U=l, V=0 in the second 
case, and finally Q=0, U=0, V=l. The value of the birefringence 
plotted in Fig. 1 is used to calculate the phase retardance at each 
wavelength A and height x, using Eq. 8. Finally, Eq. 12 is solved 
for a spectral range from 450 to 750 nm with wedges of 3-mm 
height {x varying from -1.5 to 1.5 mm). The 2D output intensi¬ 
ties, with respect to the wavelength and the height x, are shown 
in Fig. 4. Three different 2D patterns are indeed obtained for the 
3 cases, corresponding to 100% Q, 100% U, and 100% V polar¬ 
ization. 

These images are compared with real data obtained in the 
laboratory in Sect. 4 . 

In practice, when real stellar observations are obtained, 
Eq. 12 is used the other way around. We collect an image that 
corresponds to a given output intensity, and need to determine 
the entrance polarization state (I Q U V) producing the 
image. This inversion can be performed through a classical 
least-squares solving method. The parameters of the systems 
(^, 0) and the variables (x. A) have to be well known in order to 
achieve a precise measurement of I, Q, U, and V over the spec¬ 
trum. 

Fig. 5 shows the simulation of the image we would obtain 
with the polarimeter described above and a high-resolution 
(R=65000) spectrometer, for a star with an effective temperature 
of Tgff=30000 K, gravity of logy=4.0, projected rotational veloc¬ 
ity i^sinz of 50 km s“^, solar abundances, a dipolar magnetic 
field with a strength of Bpoi=1000 G inclined by j6=90° from 
the rotation axis. The star is seen equator-on (z=90°) with the 
magnetic pole facing the observer (magnetic phase is zero). The 
black lines seen in Fig. 5 are absorption lines due to the photo¬ 
sphere of the star. As the QUV Stokes parameters of the stars 
are very small compared to its intensity, no polarization modu¬ 
lation can be seen by eye. However, the information is hidden 
in the intensity, following Eqs. 8 and 12, which can be summa¬ 
rized by: 

= ^/(x. A).1(A) + xq(x, A).Q(A) + 

xlz(x, A).L[(A) + xy(x, A).y(A) (13) 


Taking the image of Fig. 5 as the observed data (correspond¬ 
ing to irt Eq. 14), we retrieve the Stokes parame¬ 

ters 1(A), Q(A), Lf(A), and V(A) of the star over the spec¬ 
trum using the least-squares method. We solve, for each wave¬ 
length, the overdetermined linear system: lout = with 
A = (xj ccq olu OLy) and S the Stokes vector we want to deter¬ 
mine (l(A) Q(A) !i(A) y(A))^. 

The results for Stokes I and Stokes V are displayed in Fig. 6. 



X 10“^ 



X (nm) 

Fig. 6. Results of the inversion code for the image presented 
in Fig. 5. Stokes I (top panel) clearly shows the photospheric 
absorption lines and Stokes V (bottom panel) shows Zeeman 
signatures at the corresponding line wavelengths, proof that the 
star hosts a magnetic field. 

The photospheric absorption lines are retrieved in the Stokes 
I spectrum. The star chosen for the simulation is magnetic, 
which means that we should find Zeeman signatures at the po¬ 
sition of the absoption lines in the Stokes V spectrum. Zeeman 
signatures are a direct proof of the presence of a magnetic field: 
the spectral lines split in wavelength in two components, which 
are left and right circularly polarized. Stokes V shows the differ¬ 
ence between these two components, and the splitting width is 
proportional to the magnetic field strength. Even if the strength 
of the Stokes V signal is 1000 times lower than the one of I, we 
retrieve the Zeeman signatures in Stokes V based on the image 
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Fig. 7. LSD profiles calculated for Stokes I and V for the syn¬ 
thetic spectrum with SNR=1000 (input, solid black line) and 
for the I and V spectra determined with the wedges concept 
and the demodulation code (output, dashed red line). They are 
strictly identical. 

simulation of Fig. 5. 

To extend the analysis further, the same simulation is per¬ 
formed but with noise in the initial synthetic spectrum of the 
star. A SNR of 1000 is used with the same star as in the previ¬ 
ous simulation. To be able to observe the Zeeman signatures, 
which have a similar amplitude as the noise, the Least-Squares 
Deconvolution (LSD) technique is applied [14]. This technique 
uses a Dirac comb corresponding to all the stellar spectral lines 
and average the Stokes I and V signals in those lines obtaining 
a mean pattern of the line profile and Zeeman signature. We ap¬ 
ply this technique for all lines from 400 to 900 nm, both directly 
to the synthetic spectrum and to the spectrum extracted from 
the 2D image. The results are plotted in Fig. 7. We retrieve a 
mean line profile and the mean Zeeman signal in Stokes V cor¬ 
responding to the input magnetic field. Moreover, the pattern 
of the LSD profiles for Stokes I and Stokes V are exactly identi¬ 
cal for the input synthetic spectrum and for the one determined 
by the demodulation code from the 2D simulation. This result 
proves the validity of the theoretical concept and of the demod¬ 
ulation code. 

D. Single component solution 

In some applications, in particular in astronomy, the amount 
of collected flux can be a crucial issue. A solution to decrease 
the number of photons lost through the polarimeter module is 
presented here. 

The previous section has shown the need for two birefrin- 
gent wedges with two different slopes to measure the full 
Stokes vector. Considering mechanical mounting issues, sym¬ 
metry around the optical axis and optical aberrations, the two 
wedges are doubled to create compact rectangular components. 
The solution proposed here is to suppress one interface for a 
total equivalent thickness. Figure 8 shows the configuration of 
the wedges. The first one, with an apex angle ^ has its fast axis 
along the horizontal direction Oy, while the second one has an 
apex angle of 2 ^ and a fast axis angle of ^ with respect to the 
Oy axis. To complete this compact single rectangular compo¬ 
nent and avoid any refraction or aberration problem, a piece of 
crystal (the same as used in the 2 wedges, MgF 2 in our case) is 



Fig. 8. Configuration of the wedges to minimize the number 
of interfaces. The first wedge has an apex angle of ^ and the 
fast axis of the crystal oriented along the Oy axis. The second 
wedge has twice that apex angle, 2 ^, with a fast axis angle of ^ 
with respect to the Oy axis. In between, we insert a piece of the 
same material (MgF 2 in our case) but with its fast axis along the 
optical axis Oz, so that is does not have any birefringent effect. 

placed in between. Its fast axis is aligned to the optical axis of 
the system Oz to suppress birefringent effects. 

Using this component, the same calculations as in the previ¬ 
ous section lead to an output intensity of: 

hut (^/ A) = 0.5 • [I + Q • cos 20 cos 20 + U • (sin O sin 20 cos 20+ 

cos O sin20) + y • (sin O sin 20 — cos O sin20 cos 20)] (14) 

with a phase retardance O of: 

0(x,A) = 2/1 tang' • . |j| ( 15 ) 

A 

Contrary to the two double-wedge system presented in 
Sect. 2A, here the value of the height x is set to 0 at the top 
of the component, where light does not cross any birefringent 
material. 

Using this single triple-component offers the same possibili¬ 
ties as the two double-wedge method detailed in Sect. 2A, but 
with the advantage of fewer components and only 4 interfaces, 
compared to 6 . This could be a precious gain in cases when the 
flux is a crucial issue, e.g., for the observation of cool stars in the 
FUV. A practical application using this component is presented 
in Sect. 3. 

3. APPLICATION TO UVMAG/ARAGO 
A. Adaptation and design 

In January 2015, the UVMag consortium [15] proposed a UV 
and visible spectropolarimetry mission, called Arago, to the 
European Space Agency (ESA). This space project offers a pre¬ 
cise polarization measurement over a high-resolution spectrum 
from 119 to 888 nm [13]. The scientific goal of Arago is to 
study the entire life of stars and planets, from their formation 
to their death, as well as feedback into the interstellar medium. 
Stellar evolution is directly driven by astrophysical processes 
such as magnetic fields and/or stellar winds. Thanks to a high- 
resolution spectrograph combined with a precised full Stokes 
polarization measurement from FUV to near-IR, Arago aims to 
obtain a full picture of the 3D dynamical environment of stars 
and planets and their interactions [16]. 
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Fig. 9. Plot of the spatial dimension x over the spectrum, for a 
given phase retardance. This corresponds to Eq. 15, with ^ = 
1.5° and O = 7r/3. As the birefringence goes to zero around 
119.5 nm, the value of x(A) plotted here goes to infinity in the 
segment [119:120] nm of the spectrum. 

Here, we present an application of the method for spectropo- 
larimetry explained in this article to the Arago instrument, 
using the single component version of the polarimeter pre¬ 
sented in Sect. 2. The system incorporates MgF 2 wedges, as 
this is the only known transparent and birefringent material 
effective over the full spectral range of Arago, from FUV to 
NIR (see Fig. 1). The apex angle ^ is set to 1.5°. Equations 14 
and 15 are used as the basis for the design of the instrument. 

We first determine the spatial period of the different modu¬ 
lation patterns that should be observed in the 2D images and 
will need to be sampled. A frequency analysis of Eq. 14 gives 
several characteristic periods of the phase retardance O. As we 
have cos O, cos 20 and cos 30 elements, the possible signals can 
have 2/1, n or 27r/3 periods of retardance. We thus need to be 
able to measure at least one period of the longest signal (27r) 
and respect the Nyquist frequency for the shortest signal (7r/3). 

The derivation of Eq. 15 for a phase retardance O = 7r/3 is 
plotted in Fig. 9. The minimum of the curve at A = 145 nm 
corresponds to the maximum of the ratio An (A) / A using the 
birefringence values presented in Fig. 1. At this wavelength, 
the space gap corresponding to a phase retardance of 7r/3 is 
65 Jim. This is the minimum value that we have to sample, i.e. 
the maximum size of the pixels. 

To find the maximum value to be sampled, i.e. the height 
of the wedges, we have to find the wavelength where the ratio 
An (A)/A is at a minimum. This occurs at the highest wave¬ 
length of study, 888 nm. The part of the spectrum between 119 
and 120 nm is not considered as the birefringence goes to zero, 
as shown in Sect. 1. At A = 888 nm, the height of the wedges, x, 
has to be large enough to observe at least one period of the low¬ 
est frequency signal, 2n. Replacing a 2/1 retardance in Eq. 15 at 
888 nm gives a minimum height of x = 2.9 mm. 

To summarize, we need to sample the height x of the wedges 
with at least 65 jim precision over more than 2.9 mm. In this 
case we do not undersample the high-frequencies in the UV and 
we observe at least one period of the low-frequency in the IR 
part of the spectrum. 

The optical design of the high-resolution spectrograph of 
Arago (spectral resolution of 35000 in the visible and 25000 in 
the UV spectral range) can be found in Pertenais et al. (2014) 
[13]. It uses a unique echelle-grating for the entire spectral 


range, and focuses all of the orders on several detectors. A 
classical CCD detector with 15 }im pixel pitch is used to image 
the 46 orders of the visible spectrum, while two micro-channel 
plate (MCP) detectors record the 178 short orders in the UV 
part of the spectrum. The pixel pitch of the MCP detectors is 
about 20 Jim. As the minimum sample we need to observe is 
65 Jim in output of the polarimeter, the spectrometer needs 
a minimum magnification factor of 1/3.25. In this way, the 
precision needed to sample correctly the signal corresponds to 
the size of a pixel. The consequence of this is that each order 
will be at least 2.9/3.25 = 0.892 mm thick. The value 900 jim is 
adopted and is the starting point to define the dimension of the 
needed detectors. For the visible detector, 46 orders of 900 jim 
with 2 black pixels (15 jim pixel pitch) between each order 
lead to a detector of at least 43 mm high (the other direction is 
defined by the main dispersion of the spectrograph). Visible 
detectors of this size are common. The UV part of the spectrum 
is spread onto two detectors, imaging 89 orders each. The 
same calculation gives two MCP detectors of at least 84 mm 
high with 20 jim pixel pitch. This is a rather large size for MCP 
detectors. 

This section has demonstrated that we can use the proposed 
method to build a static spectropolarimeter, working over a 
wide spectral range, providing all Stokes parameters at once, 
and adapted to space constraints (single static component). One 
of the drawbacks of the technique is the need for a large sensi¬ 
tive area in the UV, due to large variation of the ratio AA(A) /A 
over the spectrum. A possible solution is to modify the echelle 
grating to decrease the number of orders in the UV, but increase 
the length of the orders. In this case, the height of the detector 
needed for the UV domain decreases, and its shape becomes 
more square. 

B. Efficiency 

To compare the efficiency of the Stokes measurement with the 
concept proposed here to other existing solutions, the formal¬ 
ism developed in del Toro Iniesta et al. (2000) [17] is used. 

At a given wavelength (or given column of the image), we 
consider that the signal is modulated in n states, where n is the 
number of resolved elements or pixels. For each pixel of the col¬ 
umn, we calculate the global Mueller matrix of the polarimeter 
(see Sect. 2B). The modulation matrix for this wavelength 0(A) 
is then built with the first lines of the n Mueller matrices. Fi¬ 
nally, an optimal demodulation matrix is calculated following: 

D(A) = (o(A)^-0(A)) •0(A)’’ (16) 

Based on the coefficients of these matrices, the value of the 
extraction efficiencies of the Stokes parameters can be deter¬ 
mined for each wavelength: 

n-EDlj (17) 

where i = [I Q U V] and n is the number of pixels used for 
each wavelength. 

The calculations are performed assuming that the modula¬ 
tion is well sampled over the whole spectrum and that at least 
one full period is recorded. We find that the efficiencies are al¬ 
most constant over the spectral range with a small sinusoidal 
modulation of about 1%. 
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Analyzer angle 0 (°) 


Fig. 10. Plot of the extraction efficiencies of the Stokes parame¬ 
ters with different polarizer angles 0. 


However, the efficiencies depend directly on the angle of the 
analyzer. Figure 10 shows the variation of the efficiencies with 
the analyzer angle 6. Efficiencies are between 30% and 80% de¬ 
pending on the choice of the polarizer angle, using a dual-beam 
system. Using a simple linear polarizer decreases their values 
by a factor of 2. 

Using the alternative method with rotating plates presented 
in Sect. 2A, we obtain efficiencies varying between 40% and 
60% over the spectrum, when using a dual-beam system. There¬ 
fore, the method presented here and the polychromatric rotat¬ 
ing modulator result in comparable efficiencies. 

Finally, the choice of the polarizer angle is a free parameter 
that can be optimised depending on the science goals. For each 
instrument the angle can be tuned to reach the maximum effi¬ 
ciency for the desired Stokes parameters (see Fig. 10), e.g. 0=16° 
if one is interested in circular polarization only (efficiency of 
81% for Stokes V) or 0=55° to optimize the instrument for lin¬ 
ear polarization (efficiency of 66% for both Stokes Q and U). 

4. LABORATORY IMPLEMENTATION 

This section describes the experimental work done to verify the 
theoretical concept described above. The theoretical and simu¬ 
lated results presented in Sect. 2 are compared to actual labora¬ 
tory observations. 

A. Optical design 

To perform the laboratory tests, the first step is to control the 
entrance light source. To do so, a polychromatic halogen lamp 
is used in combination with an integrating sphere to ensure the 
use of unpolarized light. The light leaves the sphere through a 
pinhole, mimicking a point source, is then collimated and goes 
through a broadband wire grid linear polarizer to control its 
polarization state. An optional achromatic quarter waveplate, 
from the compagny Bernhard Halle Nachfolger GmbH, can be 
placed to create circular or elliptical polarization states. Finally, 
a filter wheel offers the possibilities to use monochromatic in¬ 
terferential filters or bandpass filters. The spectral region of in¬ 
terest for the laboratory validation is limited to a part of the vis¬ 
ible from about 495 to 725 nm. This spectrum, although wide, 
is shorter than the one of Arago, for example. However, these 
tests are performed only to prove the validity and demonstrate 
a practical application of the proposed concept. As explained in 
the previous section, the entire instrument can be designed to 
avoid undersampling at the edges of the spectrum, and permit 
the measurement of a much wider spectral range. 



Fig. 11. Picture of the complete optical bench showing all the 
elements, from the light source (on the right side) to the detector 
(on the left side). 


To build the polarimeter, two double wedges of MgF 2 are 
used as described in Fig. 3. The wedges, manufactured by the 
Karl Lambrecht Corporation, have submit angles of 1.5° and 
3°, a clear aperture of 10 mm, and fast axis angles of [^; 

0 ; To complete the polarimetric module of the prototype, 
a broadband wire grid polarizer is added after the wedges to 
analyse the polarization state. Its angle is chosen to be 0 = 64° 
to (1) avoid particular configurations (see end of Sect. 3B), and 
(2) obtain 3 different basic 2D patterns for Stokes Q, U and V. 
Indeed, at some positions of the polarizer, the 3 basic patterns 
can be mathematically orthogonal for example. By chosing 0 = 
64°, we consider the most general case. 

After the polarimetric module, a 50 mm focal length photo¬ 
graphic lens collimates the light coming from a slit placed right 
after the birefringent wedges and illuminates a transmissive 
grating with 600 grooves/mm to disperse the light and create 
the spectrum. As we do not have any flux problem in the lab¬ 
oratory, the slit is illuminated with a circular collimated beam, 
even though most of the light is then lost. Another solution 
would be to place the slit at the focus of a cylindrical mirror or 
lens. A second 50 mm focal length photographic lens finally 
focuses the light onto a CCD detector (6.45 }ivcv pixel pitch and 
active area of 8.67 x 6.60 mm^). The global view of the optical 
bench is shown in Fig. 11. 

B. Results and comparison to the simulations 

The experimental setup produces a 2D image. This 2D frame 
is composed of the spectrum in the dispersion direction of the 
grating (orthogonal to the height of the wedges) and the polar¬ 
ization modulation in the perpendicular direction. By rotating 
the linear polarizer placed after the light source, the frequencies 
of the intensity modulation change and we can easily recognize 
the modulation patterns of the simulations presented in Fig. 4. 

In particular, when we input 100% (Q, U, or V) polarized 
light, the specific patterns are easily recognized. However, to 
be able to measure precisely a polarization state, the instrument 
must be perfectly calibrated. 

Interferential filters with different bandwidths and central 
wavelengths are used to calibrate the position of the spectrum 
in wavelength on the detector. An example of an image ob¬ 
tained with a narrow-band filter is shown in Fig. 12. 

To calibrate the orthogonal direction of the spectrum, the 
height of the slit (3 mm) is used. In theory, the image on the 
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Fig. 12. Image obtained with a Q entrance polarization state 
and a narrow-band filter with peak-wavelength at 560 nm and 
FWHM of 10 ± 2 nm. This image is used to perform the wave¬ 
length calibration of the spectrum. 


Simulations Measurements 



Fig. 13. Comparison between the theoretical simulations (left 
panels) described in Sect. 2 and the images obtained in practice 
(right panels) for an entrance polarization of 100% Q (top), 100% 
U (middle) and 100% V (bottom).. 


detector corresponds to a variation of the height x from -1.5 to 
1.5 in Eq. 8. For processing reasons (edge effects) and vignetting 
issues, we start the observation after -1.5 mm and stop before 
1.5 mm. We thus cannot simply assume that the 2D frame goes 
from X = —1.5 to X = 1.5. Instead, the frame obtained in the lab¬ 
oratory is compared (image substraction and cross-correlation) 
with a theoretical image at given x^m and Xmax values. The 
^min and Xmax values corresponding to maximal correlation are 
extracted and the full frame is calibrated with these values. 

Once the setup is calibrated, the entrance polarization state 
is modified with the linear polarizer and the achromatic quarter 
waveplate, and the output images are recorded. The images we 
obtained are compared with the simulations and presented in 
Fig. 13. 

The demodulation code mentioned in Sect. 2 is used to mea¬ 
sure the Stokes parameters from the 2D images. For a linear en¬ 
trance polarization we are able to measure the polarization an¬ 
gle with a precision of ±1° all over the entire spectral range. A 
plot of the polarization angle value over the spectrum is shown 
in Fig. 14, for a 100% Q entrance polarization state (i.e., theori- 
cally, a 0° polarization angle). 

In a practical implementation, the final goal of such an in- 



Wavelength (nm) 


Fig. 14. The polarization angle measured in the laboratory is 
plotted for the whole spectral range for a linear polarization of 
100% Q. The measured angle is 0° ± 1°, which corresponds to a 
measurement error of 1.1%. 


strument is to measure the value of the Stokes parameters over 
the spectral range. To compensate flat-field and illumination 
issues, the actual values we determine are Q/I, U/I and V/I 
for each wavelength of interest. The plots in Fig. 15 show the 
results measured with the optical bench described above. 

The type of entrance polarization state is well retrieved us¬ 
ing the demodulation code based on the least-squares method. 
However, the amount of polarization (100%) is not properly re¬ 
trieved. 

C. Discussion 

As shown in Fig. 15, the input Stokes parameters are retrieved 
with our laboratory setup, but not perfectly. In theory we 
should measure 100% Q, U, or V, over the whole spectral range. 

First, the value determined is less than 1 (100%). This is not 
critical for polarization angle measurements defined by half of 
the tangent of the ratio between U and Q, but it is important for 
the determination of the absolute degree of polarization. This 
is mainly due to the fact that the amplitude of the measured in¬ 
tensity modulation is not maximal. Because of noise, bias, flat- 
field issues, and performance of the wedges (precise position¬ 
ing, fabrication imperfections, or non-homogeneity) the "white” 
stripes of the modulation are not equal 1 and the "black" ones 
not equal to 0. This loss of contrast propagates in the equations 
and causes this offset in the absolute Stokes parameter values. 

In addition, the determined values of the Stokes parameters 
are not constant over the spectrum. This is not due to the po¬ 
larization performances of the components, since they were all 
chosen to be achromatic. Instead, the demodulation process in¬ 
troduces some errors. The analytic solution to the linear prob¬ 
lem (Eq. 12) is used to demodulate the signal. However, some 
uncertainties apply and slightly modify this analytic solution: 
the detector is not homogeneously illuminated, the birefrin¬ 
gence curve depends on the material sample and on temper¬ 
ature, the angle of the polarizer 0 is known with 1° precision 
with respect to the entrance polarization state, the apex angle of 
the wedges ^ and the orientation of the fast axis in each wedge 
are also not perfectly known. 

While the measurements are not perfect, due to the above 
listed limitations, our experimental setup shows that the con¬ 
cept works and can produce usable results. The current proto¬ 
type will be difficult to improve because most of the limitations 
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Fig. 15. Measurements of the normalized Stokes parameters 
Q/I, U/I and V/I over the spectrum, corresponding to an en¬ 
trance light with linear polarization state of 100% Q (top), 100% 
U (middle) and 100% V (bottom). 


depend on wavelength, time, and temperature, or are inherent 
to the components. However, a more precise instrument could 
be built for space missions following this concept. 

5. CONCLUSIONS 

The use of this concept offers many advantages that can be ap¬ 
plied to various domains. As shown in this article, the idea to 
spatially modulate the polarization information makes it possi¬ 
ble to build a spectropolarimeter without any moving parts and 
which provides all Stokes parameters simultaneously. Further¬ 
more, the dependence on the wavelength is taken into account, 
and there is, theoretically, no limitation to the spectral range. 
The material used for the wedges (through its transparence and 
birefringence characteristics), in combination with the choice of 
detectors, will ultimately define the minimum and maximum 
wavelengths that can be observed. We find that MgF 2 is best 
suited for a wide range from FUV to IR. 

The optical bench we set up proves the validity of the the¬ 
oretical concept. The observational results presented here are 
comparable to the computed simulations and we retrieve the 
values of all the Stokes parameters over a given spectral range. 
The angle of linear polarization is for instance determined with 
a ±1° precision. 

This concept has thus been demonstrated and could be used 
in future space missions. 
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